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SUMMARY 


An experimental heat-transfer investigation was conducted using metallic 
tiles arranged in two staggered tile arrays with longitudinal and transverse 
gaps between the tiles. This tile arrangement results in longitudinal gaps 
terminating at intersections with transverse gaps. A direct impingement of the 
longitudinal gap flow on the forward face of a downstream blocking tile results; 
hence, high localized heating occurs at this impingement region. The tile arrays 
were tested (1) to define the impingement heating distributed on^a downstream 
blocking tile and (2) to determine the influence of tile and gap geometries in 
sufficient detail to develop empirical relationships for impingement heating in 
laminar and turbulent boundary layers. The tests were conducted in the Langley 
8-foot high-temperature structures tunnel at a nominal Mach number of 7, a nom- 
inal total temperature of 1800 K, and free-stream unit Reynolds numbers from 
1.0 X 10® to 4.8 X 10® per m. Tile forward-facing wall slopes of 60°, 75°, and 
90° were tested along with longitudinal gap lengths of 15.24 and 30.48 cm and 
gap widths of 0.10, 0.18, 0.30, and 0.41 cm. The boundary-layer displacement 
thickness covered a range from 0.36 to 1.62 cm. No laminar data were obtained 
on the 90° tile array. Flat-plate calibration tests provided comparison data. 

For laminar boundary layers, the heating rates observed on the downstream 
blocking tile in the transverse gap direction for gap widths of 0.10 and 0.18 cm 
indicate that the gap flow is primarily two-dimensional in nature (i.e., flow in 
the longitudinal gap does not turn and flow in a transverse gap). Thus, the 
longitudinal gap flow has a negligible effect on the heating rates on the block- 
ing tile for gap widths less than 0.18 cm.. For a gap width of 0.41 cm, the 
heating rates in the transverse direction are less uniform. This lack of uni- 
formity suggests that the basic flow is three-dimensional and that the flow in 
the longitudinal gap is becoming more influential. The data indicate that the 
maximum impingement heating of 2.4 times the flat-plate heating rate occurs for 
a gap width of 0.41 cm, and for the two smaller gap widths, the impingement heat- 
ing was less than the surface heating. For turbulent boundary layers, the flow 
in the transverse gap is primarily three-dimensional in nature. The region 
affected by the three-dimensional flow extends at least three gap widths in the 
transverse and longitudinal directions. The maximum hearing rate observed was 
4.5 times the surface heating rate and occurred for a forward-facing wall slope 
of 90° and a gap width of 0.41 cm. Changing the slope of the forward-facing 
wall from 90° to 60° decreased the heating substantially in the impingement 
region. Decreasing the gap width also causes a decrease in impingement heating. 
In both laminar and turbulent flow the presence of gaps influences surface 
heating . 

Empirical relationships to correlate Impingement heating with gap flow and 
geometric parameters for laminar and turbulent flow were developed through com- 
bined regression analysis and graphic techniques. These relationships accu- 
rately predict the impingement heating rate and the influence of gap flow and 
geometric parameters within the data range used to develop the relationships. 



INTRODUCTION 


The thermal protection system (TPS) of the space shuttle orbiter consists 
of a silica-based material for the reusable surface insulation (RSI) . The RSI 
is bonded to a strain isolator pad which is bonded to the primary structure and 
must withstand temperatures of approximately 1700 K on the bottom surface of 
the orbiter and maintain a primary structure temperature below 450 K during 
entry. Arranged in a staggered pattern, the RSI tiles are applied to the sur- 
face with gaps between tiles to accommodate thermal expansion and contraction 
and mechanical deflections of the underlying structure as well as to allow for 
thermal expansion of the tile material. The tiles are nominally 15 by 15 cm 
and the thickness (0.50 to 9.00 cm) is varied according to the expected heat 
load. The gaps can locally disrupt the external boundary layer and are, there- 
fore, sources for high localized heating during entry into the Earth's 
atmosphere . 

Past experimental studies (refs. 1 to 3) have identified as a potential 
problem high localized heating regions on RSI tiles where the flow in a longi- 
tudinal gap intersects a transverse gap and impinges on a downstream blocking 
tile. The intensity of the localized heating is influenced by tile and gap 
geometries, local flow conditions, and tile orientation to the direction of 
flow. However, these studies did not address the localized heating regions in 
sufficient detail to permit accurate assessment of the overall impact of this 
impingement heating at the downstream blocking tile. Therefore, an experimen- 
tal heat-transfer investigation was conducted to define the heating distribu- 
tion in the impingement region and to determine the influence of tile and gap 
geometries in sufficient detail to develop empirical relationships for impinge- 
ment heating in laminar and turbulent boundary layers. 

This study was carried out using metallic tiles in staggered arrays. The 
test parameters were longitudinal gap length, gap width, slope of the forward- 
facing walls, boundary-layer thickness and state (laminar/turbulent) , and 
Reynolds number. The tiles were densely instrumented in areas of expected high 
heating to obtain well-defined impingement heating distributions. The cold-wall 
heating rates were obtained in the Langley 8-foot high-temperature structures 
tunnel at a nominal Mach number of 7, a ncaninal total temperature of 1800 K, 
free-stream unit Reynolds numbers from 1.0 x 10® to 4.8 x 10® per m, and a free- 
stream dynamic pressure range of 15 to 61 kPa. 


SYMBOLS 

The units used for the physical quantities in this paper are given in the 
International System of Units (SI) . U.S. Customary Units were used for the 
principal measurements and calculations. 

Cp specific heat, J/kg-K 

d total tile thickness, cm 

dT/dt temperature rise rate, K/s 
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L length of longitudinal gap (fig. 3), cm 

I total length of calibration panel (fig. 11 )f cm 

p^^^ total pressure in combustor, MPa 

p pressure, MPa 

q dynamic pressure, kPa 

q cold-wall heating rate, kW/m^ 

R unit Reynolds number, per m 

R^ local Reynolds number based on distance from panel holder leading 

edge 


s 

l^t, c 
w 


local Reynolds number based on gap length 
surface distance from tile radius mid arc, cm 
total temperature in combustor, K 
gap width (fig. 3), cm 

distance from leading edges of calibration panel (fig. 11), cm 
vertical distance into gap (fig. 2) , cm 
angle of attack, deg 

boundary-layer displacement thickness, cm 
slope of forward-facing wall (fig. 3), deg 
slope of forward-facing wall, rad 
density, kg/m^ 
material thickness, cm 
theoretical heating-rate ratio 
Subscripts: 
atm atmospheric 


a 

6 * 

0 

0rad 

P 

T 

4 > 


em 


FP 


experimental 

empirical 

flat-plate condition 



I 


impingement 


S surface 


TEST APPARATUS 
Model Description 

Calibration panel *- As shown in figure 1, the calibration panel was instru- 
mented with 35 number 30-gage chrome 1-alumel thermocouples which were spot 
welded to the back surface of a 0, 05-cm-thick 304 stainless-steel plate. (See 
ref. 1 for further details.) The calibration panel was used to determine tunnel 
operating parameters which would yield a laminar boundary layer over the test 
panels and to obtain the corresponding cold-wall laminar heating rates for ref- 
erence data. 

Tile array panels .- The tile arrays shown in figure 2 were nominally 46 by 
48 by 6,35 cm and consisted of solid aluminum tiles surrounding two instrumented 
stainless-steel thin-wall (0.08 cm) tiles. Metallic tiles were used to simulate 
the RSI tiles because they were easier to handle, and the tile material has no 
effect on the cold-wall heating rates. Each tile had an edge radius of 0.25 cm. 
The model schematics indicate orientation to the flow direction, instrumented 
tile locations, and geometric parameters. The geometric parameters listed are 
longitudinal gap length L, gap width W, and slope of the forward-facing wall 
9, The test ranges of these parameters are compared with the nominal shuttle 
geometric parameters in figure 2(c). The forward-facing walls of the thin- 
walled tiles are sloped at 60^, 75^, or 90^. For convenience, the 60^ and 75® 
walls were provided on the same tile (figs. 2(a) and 2(b)) and required only 
180® rotation and appropriate filler bars to change configurations. All the 
tiles were mounted on a 0. 32-cm-thick stainless-steel plate with slotted holes 
to permit adjustment of tile position to change gap width. The gap widths were 
0.10, 0.18, 0.30, and 0.41 cm. liOngitudinal gap lengths were 15.24 and 30.48 cm. 
photographs of the models are shown in figure 3, 

The thin-wall tiles were instrumented with 30-gage chromel-alumel thermo- 
couples spot welded to the inside surface of the tiles. The thermocouples were 
distributed both longitudinally and laterally over the tile surface as shown in 
figures 4 and 5. Thermocouples are concentrated in the impingement regions at 
the downstream end of the longitudinal gap to define the expected highly nonuni- 
form heating distributions. Each 60®/75® tile had 43 thermocouples and each 90® 
tile had 25 thermocouples. The numbers in the figures indicate thermocouple 
designations for tiles at the end of the 15.24-cm longitudinal gap. Thermocou- 
ple designations for L = 30,48 cm may be obtained by adding or subtracting the 
constants indicated in the respective figures. 

In order to define heating distributions in the impingement region, pre- 
cise thermocouple locations had to be known. The thermocouple locations were 
precisely determined by using optical measurement equipment and radiographs of 
the thin-wall tiles to identify the mean longitudinal center line of the thermo- 
couples. These thermocouples varied within ±0,03 cm of the scribed mean longi- 
tudinal center line. A precision alignment rig was used to align the mean ther- 
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inocouple center line with the upstream longitudinal gap center line and to set 
the desired gap widths. 


Panel Holder 

The models were mounted in a panel holder (figs. 6 and 7) which can accOTi- 
modate test panels up to 152 by 108 cm. (See refs. 4 and 5.) The aerodynamic 
surface is covered with 1 .25-cm-thick low conductivity Glasrock tiles. Aerody- 
namic fences (fig. 6) provide uniform two-dimensional flow over the entire aero- 
dynamic surface. A blunt leading edge with a radius of 0.95 cm is used on the 
panel holder to promote a laminar boundary layer ^ and a sharp leading edge with 
a lateral row of boundary-layer trips is used to prc^ote a turbulent boundary 
layer over the aerodynamic surface of the panel holder. Two test locations 
were provided in the panel holder 117 and 188 cm from the leading edge. The 
calibration panel and the 60^/75® sloping tile array were located in the forward 
location (location I) , and the 90^ tile array was located in the rear location 
(location II) . 


Facility 

The tests were performed in the Langley 8-foot high-temperature structures 
tunnel (fig. 8) which is a large blowdown facility that simulates aerodynamic 
heating and pressure loading at a nominal Mach number of 7 and an altitude 
between 25 and 40 km. The high energy needed for this simulation is obtained 
by burning a mixture of methane and air under pressure in the combustor and 
expanding the products of combustion through a conical-contoured nozzle into 
the open jet test chamber. The flow enters a supersonic diffuser where an air 
ejector pumps the flow through a mixing tube and exhausts the flow to the atmo- 
sphere through a subsonic diffuser. This tunnel operates at a combustor total 
temperature T^.^^ between 1400 and 2000 K and ccMnbustor total pressures Pt,c 

from 4.1 to 24.1 MPa. Free-stream unit Reynolds numbers range from 1.0 x lo^ to 
10.0 X 10^ per m. 

Test models are mounted on an elevator and inserted into the test stream 
after the test conditions are established to avoid exposing test models to tun- 
nel transients which occur during tunnel start-up and shutdown. A model pitch 
system provides an angle-of-attack range of ±20^. More detailed information 
can be found in references 4 and 5. 


TESTS 

Wind-tunnel tests were made on a calibration panel and on tile array pan- 
els. The calibration panel was used to determine tunnel flow conditions which 
would yield a laminar boundary layer over the test panels. The calibration 
panel was also used to obtain corresponding cold-wall laminar heating rates for 
data correlation with nonsmooth surfaces. The tile array panels were used to 
define the heating distribution in the impingement region and to determine the 
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influence of tile and gap geometries in sufficient detail to develop empirical 
relationships for impingement heating in laminar and turbulent boundary layers. 


Test Procedures and Data Reduction 

The test procedures were to establish the tunnel equilibrium flow condi- 
tions and then insert the model into the test stream. The model was pitched 
during insertion to the desired angle of attack to give it, as nearly as possi- 
ble, a step-function exposure to the test condition. A typical insertion his- 
tory showir^ the model position, the angle of attack, and a typical temperature 
response on the metallic tiles is shown in figure 9, The model traverses 210 cm 
to the tunnel center line in a little over 2 s. The panel holder reaches the 
desired angle of attack at approximately 1 s after insertion. The tile tempera- 
ture response shows that the maximum temperature rise rate occurs before the 
model reaches the tunnel center line. However, until the model reaches the tun- 
nel center line, the entire panel holder is not exposed to a uniform test envi- 
ronment. Therefore, the time of data analysis, as shown in figure 9, was chosen 
as the time when the model reached the tunnel center line. 

Models and tunnel instrumentation data were recorded by high-speed digital 
recorders at 20 and 200 frames per s. Cold-wall heating rates were calculated 
from these outputs using the one -dimensional transient heat balance equation 
(q = pcpT (dT/dt>) . The temperature rise rate dT/dt was calculated by averag- 
ing the model .temperatures over 0.25— s intervals and then determining the temper- 
ature rise rate between each interval. Therefore, the heating rates were deter- 
mined each 0.25 s. 

The one-dimensional heat balance equation equates the convective energy 
entering the surface to the energy stored and does not consider the effect of 
conduction and radiation. These assumptions are considered reasonable since 
the temperature— time slopes were taken early in the tests when the surface tern- 
pgratures were relatively low. In addition, calculations made to account for 
conduction and radiation effects indicated these heat losses were less than 
2 percent at the time data were taken. 


Conditions 


The tunnel conditions for the calibration series are given in table I. 

The calibration panel was tested three times in ncxninal tunnel conditions of 
= 1800 K, dynamic pressures of 14 and 42 kPa, free-stream unit Reynolds 
numbers of 1 .0 x lO^ and 3.0 x lO® per m, and a = 0° and 4.5°. The three 
selected combinations of these tunnel conditions resulted in laminar flow over 
the calibration panel. 

Tunnel conditions and test parameters for the tile array models are given 
in table II. The tile arrays were tested a total of 17 times; 7 tests were 
made in a laminar boundary layer and 10 in a turbulent boundary layer. The 
laminar tests were at p = 1800 K (nominal), dynamic pressures of 15 and 
32 kPa, free-stream unit Reynolds numbers of 1.0 x lO and 2.6 x io per m, and 
Qj - 0°. The turbulent boundary-layer tests were at ft,c ~ 1800 K (nominal). 


6 


dynamic pressures of 25 and 61 kPa, free-streara unit Reynolds numbers of 
2.0 X 10® and 4.0 x 10® per m, and a = 0° and 7.5®. These conditions yielded 
laminar 6* values of 0,36 and 0.57 cm and turbulent 6* values between 0.85 
and 1.62 cm (ref. 6). The free-stream tunnel conditions were determined frcan 
temperatures and pressures measured in the ccwnbustor and are based on the 
thermal, transport, and flow properties of methane-air combustion products as 
reported in reference 7 and in the tunnel surveys of reference 5. 


RESULTS AND DISCUSSION 
Flow Conditions and Test Data 

Calibration test data .- Typical experimental cold-wall convective heating 
rate distributions over the calibration panels are compared with theory in fig- 
ure 10 for laminar and turbulent boundary layers. The laminar data were 
obtained on the calibration panel in the forward region of the panel holder. 

(See inserts in fig. 10(ay.)^ The turbulent data were obtained from an earlier 
study using a large calibration panel (ref. 5) which covered the entire test 
area of the panel holder. A sharp leading edge with a lateral row of spherical 
boundary-layer trips (see inserts in fig. 10(b)) was used to promote turbulence. 
The theoretical predictions shown by the dashed curves were calculated using 
equations and plots from reference 5. In general, the laminar boundary-layer 
theory was 10 percent higher than the laminar experimental data (fig. 10(a)) 
obtained on the small calibration panel. Reference 5 reported turbulent 
boundary-layer theory to be approximately 30 percent higher than the interpo- 
lated experimental data (fig. 10(b)) which was an expected error. The inter- 
polated experimental data were obtained by interpolating between data in 
reference 5. 

Laminar flow conditions existed over location I for the range of test 
conditions listed in table I. Flat-plate heating rates in laminar flow range 
between 6.5 and 9.0kW/m2. Results from reference 5 showed that uniform turbu- 
lent flow could be maintained over the entire test surface of the panel holder 
for test conditions of 6.8 MPa ^ Pt,c ^ 17.3 MPa and 0® S a $ 15®. The flat- 
plate heating rates in turbulent flow range between 63 and 137 kw/m^. 

Tile array test data .- Experimental tile heating rates nondimensionalized 
by flat-plate heating rates are tabulated in table III for each thermocouple 
location. Flat-plate heating rates qpp in a laminar boundary layer were 
obtained on the small calibration panel at the location of the thin-wall tile 
centers except for the condition Pt,c “ 10.2 MPa and a = 0.3®. At this test 
condition the theoretical value of the qpp was used because calibration data 
were not available. This theoretical value is considered to be reliable because 
of the good agreement between experiment and theory shown in figure 10(a). As 
previously mentioned, the flat-plate heating rates in a turbulent boundary layer 
were obtained from the experimental data of reference 5 at the location of the 
thin-wall tile centers. 


Under certain flow conditions the presence of the tile array (location I) 
appeared to create transitional flow where laminar flow existed over the cali- 
bration panel. For example, whereas laminar heating rates were obtained on the 
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smooth calibration panel at ^ = 4,0 MPa and oc = 4.5® (see table 1)'^ heat- 
ing rates indicative of transition were obtained on the tile array at the same 
test conditions and location. This phenomenon was also noted in reference 8. 
Heating rate data for location II were indicative of transitional flow at all 
flow conditions with the arrays present. (In general for the flow conditions 
presented, ref. 5 indicated laminar flow for a smooth surface.)* No laminar 
heating rate data were obtained for the tile array in location II. Tunnel con- 
ditions, test parameters (table II), and heating rate ratios (table III) are 
tabulated only for runs where the boundary layer was either laminar or turbulent. 


Aerodynamic Heating in Laminar Flow 

Tests were conducted under laminar flow conditions to identify local flow 
characteristics and to determine the magnitude and distributions of the convec- 
tive heating in the impingement region downstream of a longitudinal gap as 
affected by changes in gap geometries. All the laminar data plotted are for 
test conditions of Pt,c = ® 

75®. 


Heating rate distributions in impingement region .- Typical laminar convec- 
tive heating rate distributions along the longitudinal center line of the thin- 
wall tiles (L = 15.24 cm) are shown in figure 11 for W = 0.18 and 0.41 cm. 

The local heating rate is nondimens ionalized by the flat-plate heating rate. 

The forward-facing walls of the tiles have been rotated up in the sketches for 
clarity in the plots. Thermocouple numbers in the sketches correspond to those 
in figure 4. The heating rate distributions in figure 11 show that the impinge- 
ment region heating rates increase with increasing W and 0 . Peak heating at 
all gap widths occurs on the top surface near the edge radius of the tiles, and 
the maximum heating rate ratio is 2.2. Results for gap widths of 0.10 and 
0.30 cm follow the same trends. 

The heating rate distributions down the forward face for all tile array 
configurations with L = 15.24 cm are shown in figure 12. Locations of the 
thermocouples are shown in the inserts. In all cases the heating decreases 
significantly with increasing depth into the gap. In fact, the maximum heat- 
ing rate was only O.lcjpp at 30 percent of the tile depth. Data shown for 
q/qpp < 0.02 are of questionable accuracy because of limitations in the data 
recording system. By integrating the center-line heating distribution on the 
forward-facing tile wall, the heat load per unit width is shown to decrease 
with decreasing 0 for W = 0.18 cm even though the surface area increases. 

The percent decrease in the impingement heating rate and in the heat load per 
unit width obtained by decreasing 0 from 75® to 60® is 31 percent and 56 per- 
cent, respectively, for W = 0.18 cm. Very little decrease is shown for 
W » 0.41 cm. 

Details of the heating in the impingement regions are shown in figure 13 
for forward-facing wall slopes of 60® and 75® with gap widths of 0.18 and 0.41 cm 
(longitudinal gap length is 15.24 cm). The short-dashed lines represent outlines 
of the upstream gap, and the long-dashed lines represent tangencies between the 
edge radius and top and forward-facing tile surfaces. Here the basic flow phe- 
nomenon created by longitudinal and transverse gaps is characterized by isomet- 
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ric heating rate contours. As these contours were generated, the heating rates 
were assumed to be symmetric about the mean thermocouple center line. 


The uniformity of the heating rates in the transverse direction for 
W = 0.18 cm indicates that the gap flow is primarily two-dimensional in nature 
and that the longitudinal gap flow has a negligible effect on these heating 
rates for gap widths less than 0.18 cm. For a gap width of 0.41 cm, the con- 
tours are less uniform, which suggests that the basic flow is three-dimensional 
and the flow in the longitudinal gap is beccMning more influential. As the gap 
width increases, the flow in the longitudinal gap also influences the heating 
on the top surfaces. Results for gap widths of 0.10 and 0.30 cm follow the 
same trends. 

Effect of tile and gap geometry .- The effects of gap width and longitudinal 
gap length on surface and impingement heating rates are shown in figure 14 (qg 
obtained from thermocouple 55; qj obtained from thermocouple 49 or 61). Data 

are shown for B = 60° and 75° and the two values of gap length (L = 15.24 and 

30.48 cm>. The general trend is an increase in surface heating with an increase 
in W, 0, and L. This increase in surface heating caused by the presence of 
gaps was also noted in reference 8. The results shown in figure 13 suggest that 
the increased heating results from changes in flow down the longitudinal gap 
due to increased gap widths. The maximum increase in surface heating of 0.5qpp 
occurred at the maximum gap width (W = 0.41 cm). However, for the smaller gap 
widths the increase in surface heating is less than 10 percent. The impingement 
heating rate qj is influenced similarly by wider gaps and longer longitudinal 
gaps. The maximum impingement heating of 2.4qpp also occurs for W = 0.41 cm. 
For the two smaller gap widths, qj was less than the surface heating rate. 

The maximum increase in qj due to gap length for 0 = 75° was 60 percent and 
occurred for W = 0.30 cm; data for 0 = 60° indicate similar trends. As pre- 
viously stated, decreasing the slope of the forward-facing wall frcm 75° to 60° 

generally decreases impingement heating and surface heating levels. 


Aerodynamic Heating in Turbulent Flow 

Past investigators of gap heating in a turbulent boundary layer have 
studied the effects of W, L, 6*, and tile orientation. However, the effect 
of sloping the forward-facing wall was not investigated nor has the impingement 
heating region been well defined. Therefore, these effects as well as the 
effects of W and L were investigated for the 60°, 75°, and 90° tile arrays. 
Test conditions for the turbulent data plotted were Pt,c “ MPa, a = 7.5°, 
and 6* = 1.01, 1.36, and 1.62 cm. 

Heating rate distributions in impingement region .- Typical nondimensional- 
ized heating rate distributions along the longitudinal center line of the tiles 
are Shown in figure 15 for forward-facing wall slopes of 60°, 75°, and 90°. 

These heating rate distributions show that the heating on the forward-facing 
wall downstream of a longitudinal gap increases rapidly to a maximum near the 
edge radius, and the heating rates decrease rapidly aft of the edge radius until 
they reach a near uniform flat-plate value at the center of the tile. For the 
indicated test conditions, peak heating on the 60° and 75° forward-facing walls 
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is approximately 2.5 to 3,0 times the flat-plate value, whereas on the 90^ 
walls, the heating is 4.3 times greater. However, this result depends on many 
other parameters discussed in later sections. 

The heating distributions down the forward-facing wall for all configura- 
tions of the tile arrays with L = 15.24 cm are shown in figure 16# All distri- 
butions show the same general trend, a constant heating rate region at the top 
which decreases rapidly with gap depth. This constant heating level occurs 
over the upper 15 percent of the gap depth for 0 = 90® with W = 0.30 and 
0.41 cm and 0 = 75® with W = 0.41 cm. For all other gap geometries the 
constant heating only occurs over the upper 5 percent of the gap depth. Both 
the maximum heating rate and the constant heating rate region increase with 
increases in W and 0 . These trends are scmiewhat different from the laminar 
trends (fig. 12) in that the laminar results did not indicate a constant heating 
rate region. At all test conditions the turbulent heating has reached O.lqpp or 
less at 60 percent of the tile depth as compared with laminar heating which 
reached 0.1 qpp at 30 percent of the tile depth. As in the laminar case, inte- 
grated heat load per unit width on the forward-facing tile wall decreases with 
decreasing 0 even though the surface area increases. The percent decrease in 
the heat load per unit width obtained by reducing 0 from 90® to 60® and 90® to 

75® with W = 0.18 cm is 23 percent and 10 percent, respectively. 

The point of peak heating shifted slightly with changing gap width, as 
illustrated by the heat rate vector plots shown in figure 17. The data shown 
are for the 90® tile array; however, the other tile array (0 = 60® and 75®) 
showed similar trends. For W * 0.10 cm the maximum heating occurred at the 
tangency point of the edge radius and the top surface. As W was increased, 
the maximum heating rate moved around the edge radius to approximately mid arc 
for W = 0.18 cm and to the tangency of the forward-facing wall and edge radius 

for W = 0.30 and 0.41 cm. 

Impingement region isometric heating rate contour plots for 0 = 60®, 75®, 
and 90®, at W = 0.18 and 0.41 cm, are shown in figure 18. The heating rates 
used in generating the contours were assumed to be symmetric about the mean 
thermocouple center line. As in figure 13, the dashed lines represent outlines 
of the upstream gap and the tangencies between the edge radius and top and 
forward-facing tile surfaces. The outlines are shown for reference only. The 
isometric heating rate contours indicate that the flow is three-dimensional. As 
was found for laminar flow, increases in W and 0 cause an increase in heat- 
ing intensity. The region affected by the three-dimensional flow extends at 
least three gap widths in the transverse and longitudinal directions. As noted 
previously, peak heating regions move down the tile face as gap width is 
increased. Results for gap widths of 0.10 and 0.30 cm follow the same trends. 
Peak heating rate level for all cases is greater than qppr and a maximum heat- 
ing rate ratio of 4.5 occurs for 0 - 90® and W = 0.41 cm. 

Effect of tile and gap geometry .- The influence of the slope of the forward- 
facing wall on impingement heating is shown in figure 19 where the impingement 
heating rate (qj obtained from thermocouple 49, 61 , or 93) nondimensionalized 
by the flat-plate heating rate value is plotted against the forward-facing wall 
slope 0 . Data are plotted for 6 * * 1.01 cm at 0 = 60® or 75®, and for 
6 * = 1.36 cm at 0 = 90®. (No data were taken on the 90® tile array for 
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6* = 1.01 cm.) Experimental data (ref. 1) indicate that heating decreases with 
increases in 6*; consequently, heating rates for 9 = 90° at 6* = 1.01 cm 
would be greater than heating rates shown for 6* = 1.36 cm. The 60° and 75° 
data indicate that 0 significantly influences impingement heating. There is 
insufficient data to show details of the trend from 60° to 90°; however, the 
influence of 9 can be considered analytically and is discussed in a later 
section. 

The effects of gap width and longitudinal gap length on impingement heating 
rate are shown in figure 20. As noted previously, the impingement heating rate 
increases with gap width and the slope of the forward— facing tile wall; the 
effect of longitudinal gap length is not as obvious. For the longer length gap 
and 0 = 60° (open squares, fig. 20(a)), increasing the gap width causes the 
impingement heating rate first to decrease and then to increase (at W > 0.18 cm) 
to levels that can be greater than the corresponding heating rate level for the 
shorter gap (open circles). The data for 9 = 75° (solid symbols, fig. 20(a)) 
appear to substantiate the trends noted above for 9 = 60°; however, for 9 = 90° 
the impingement heating rate for the longer length gap does not increase to 
levels higher than the corresponding levels for the shorter length gap. This 
inconsistency with data at 9 = 60° and 75° is not fully understood because of 
limited data; however, the inconsistency may be attributed to the thicker bound- 
ary layer v^ich is the only difference in the data. Unpublished results from 
tests at NASA Langley and Ames Research Centers indicate that a periodic heat- 
ing pattern exists along the length of longitudinal gaps (gap lengths were 15, 

48, and 84 cm with gap widths up to 0.48 cm). The fluctuations in the impinge- 
ment heating rate indicated in figure 20 may be due to the same unexplained 
phenomenon noted in the longitudinal gaps. 

Geometric effects of 0, W, and L on the midpoint surface heating rates 
are indicated in figure 21 . The geometric parameters have an overall opposite 
effect on surface heating than on impingement heating; that is, if a change in 
a parameter tends to increase the impingement heating, then that change tends 
to decrease the surface heating. However, the same periodic heating rate effect 
is evident, except that the shorter longitudinal gaps also exhibit the effect. 

It should also be noted that the periodic surface heating rate is in phase with 
the periodic impingement heating (e.g., a decrease in impingement heating is 
accompanied by a corresponding decrease in surface heating) . 


Correlation of Laminar and Turbulent Data 


Past investigators (ref. 9) have proposed a correlation parameter which 
relates heating in the impingement region with boundary-layer and tile gap geo- 
metric parameters. The authors of reference 10 extended and refined the corre- 
lation to apply to a wider range of test data. The intent was to define an 
impingement heating envelope for the maximum heating at any arbitrary flow and 
gap conditions. The correlation parameter [|(6*/W) 2 (z/L)^ relates the impinge- 
ment heating rate to boundary-layer displacement thickness 6*, gap width W, 
data point vertical location z, and longitudinal gap length L. This correla- 
tion parameter does not adequately correlate the laminar data. Therefore, the 
parameter was modified to include laminar data by including a ratio <]> where 
4> is the ratio of the theoretical heating rate for turbulent flow divided by 
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the theoretical heating rate for the actual flow conditions over a smooth sur- 
face* (The ratio equals one when the actual flow is turbulent and is greater 
than one when the actual flow is laminar*) The present data and turbulent flow 
data from references 1 and 9 plotted against the modified correlation parameter 
are compared with the upper bound for impingement heating from reference 10 in 
figure 22* The correlation parameter gives only fair agreement with experimen- 
tal data and does not identify all the governing parameters and their degree of 
influence* Consequently, empirical relationships were developed which could 
accurately predict the maximum impingement heating as well as indicate the gov- 
erning parameters and their degree of influence* 

The experimental data presented herein indicate that impingement heating 
is dependent on governing parameters of gap geometry (W, L, and 0) and flow 
conditions (boundary-layer thickness and Reynolds number) . Even though data 
from references 1 and 9 were included in the data base to enlarge the parameter 
range, sufficient data were not available to isolate the independent effects of 
all the parameters* This deficiency can be bridged by using a regression anal- 
ysis* The regression analysis method of references 11 and 12 was used to inves- 
tigate many functional relationships between the heating rate data and the flow 
and geometric parameters* 

An equation of the general form 


y = ao xf xf3 . . . ( 1 ) 

was found to yield the "best fit" to the data. By studying the experimental 
data, the physical problem, and the computed statistical information about each 
parameter and functional relationship, parameters having a significant influence 
on the impingement heating rates could be determined. The following table shows 
all the parameters investigated and the parameters that were found to exert the 
most significant influence on the impingement heating rates; 


Boundar y-1 aye r 
state 

Parameters investigated 

Parameters of 
significant influence 

Laminar 

6*, 6*2, w, w2, 6*/W, (6*/W)2, 

z, W/6*z, (W/6*z)2, (W/6*z)^, 

(W/6*z)0*5, L, z/L, (|), In (|>, 
^i,Lf sin 0, Grad 

6*/W, Rz/ sin 0 

Turbulent 

6*, 6*2, w, w2, 6*/W, (6*/W)2, 

z, W/6*z, (W/6*z)2, (w/6*z)^, 

(W/6*z)0*5, L, z/L, <|), In 4), 
^L' sin 0, ^rad' 

Ps/Patm 

Ps/Patm» ^Z,L' ® 


For a laminar boundary-layer state, all the expected governing flow and geomet- 
ric parameters had a significant influence. Some of the parameters are interre- 
lated (e.g., and indicating an interaction between the flow and 
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geometric parameters. The functional relationship of 0 could not be adequately 
defined as an independent parameter. Consequently, the effect of 6 was quali- 
tatively taken into account by sin 0 and included with the dependent variable 
as q/qpp sin 0. For a turbulent boundary— layer state the same governing flow 
and geometric parameters were found to be influential except for Rj. As in 
the laminar case, Rj was expected to be an influential parameter; however, 
turbulent data correlated better using the surface pressure pg. Although there 
is a direct relationship between pg and Rj, the latter parameter also reflects 
changes in other parameters (e.g., velocity for the present tests); in addition 
the turbulent gap flow phenomenon is apparently better characterized by the sur- 
face pressure pg. 


In obtaining the regression equation, all data were initially input into 
the regression analysis, and those observations with high residuals were then 
eliminated to improve the overall fit (for laminar flow, 10 of 14 observations 
were retained; for turbulent flow, 27 of 51 observations were retained) . The 
final data used in obtaining the regression equations are listed in table IV. 
The observation number for each of the first 30 observations indicates the run 
number (first two digits! and the thermocouple number (last three digits)' for 
the present data. The resulting regression equations are as follows: 


Laminar : 


q 

qpp sin 0 


/6*\-l .12 
= 17 000 — 

\W / 


R 


-0.07 


r-0.53 


(2) 


Turbulent; 


q 

qpp sin 0 



6*\-0.44 


p^+0 • 05 

I fit 


^ Ps 

^ I 

Patm;^ 


0.18 


(3) 


These equations were then compared with the test data to assess their adequacy 
in predicting the impingement heating level and the influence of the governing 
parameters. The equation for turbulent flow did not fully account for the 
effects of gap width and slope of the forward-facing wall. Therefore, the 
regression equation (eq. (3)) was used as the basis to graphically isolate the 
effects of W and 0 . This combined regression analysis and graphic technique 
resulted in the following equation; 



The form of the additional terms in equation (4) , which resulted from the 
graphic technique, was inconsistent with the form requirement of the regression 
analysis technique. (See eq. (1).) The adequacy of the combined equation 
(eq. (4)) is illustrated in figure 23. The regression equation for turbulent 
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flow (eq. (3)) does not predict the heating levels for the effects of 9 and 
W (fig. 23(a)); however, the combined equation (eq. (4)) accurately predicts 
both the trend and levels. The data indicate a reduction in heating of approx- 
imately 32 percent when reducing 9 from 90® to 60®; however, equation (4) 
indicates that a reduction in heating of 38 percent may be possible. The com- 
bined equation also improved predictions for the effect of W as shown in 
figure 24(b). 

The accuracy of these empirical relationships (eqs. (2) and (4)) is illus- 
trated in figure 24. Ninety percent of the laminar and turbulent impingement 
heating rates can be predicted within ±29 percent. Consequently, equations (2) 
and (4) accurately predict the impingement heating rate and the influence of 
gap flow and gecxnetric parameters within the data range used to develop the 
relationships. 


CONCLUDING RQ!ARKS 

An experimental heat transfer investigation was conducted using metallic 
tiles arranged in two staggered tile arrays with longitudinal and transverse 
gaps between the tiles. This tile arrangement results in longitudinal gaps ter- 
minating at intersections with transverse gaps \rfiich creates a direct impinge- 
ment of the longitudinal gap flow on the forward face of a downstream blocking 
tile; hence, high localized heating occurs at this impingement region. The tile 
arrays were tested to define the impingement heating distribution on a down- 
stream blocking tile and to determine the influence of tile and gap geometries 
in sufficient detail to develop empirical relationships for impingement heating 
in laminar and turbulent boundary layers. The tests were conducted in the 
Langley 8-foot high-temperature structures tunnel at a nominal Mach number of 7, 
a nominal total temperature of 1800 K, and free— stream unit Reynolds numbers 
from 1.0 X 10® to 4.8 x lO® per m. Tile forward-facing wall slopes of 60®, 75®, 
and 90® were tested along with longitudinal gap lengths of 15.24 and 30.48 cm 
and gap widths of 0.10, 0.18, 0.30, and 0.41 cm. The boundary-layer displace- 
ment thickness covered a range from 0.36 to 1.62 cm. No laminar data were 
obtained on the 90° tile array. 

For laminar boundary layers, the heating rates observed on the downstream 
blocking tile in the transverse gap direction for a gap width of 0.18 cm indi- 
cate that the gap flow is primarily two-dimensional in nature (i.e., flow in 
the longitudinal gap does not turn and flow in a transverse gap) . Thus, the 
longitudinal gap flow has a negligible effect on the heating rates on the block- 
ing tile for gap widths less than 0.18 cm. For a gap width of 0.41 cm the iso- 
metric heating rate contours are less uniform and suggest that the basic flow 
is three-dimensional; thus, the flow in the longitudinal gap is becoming more 
influential. The data indicate that the maximum impingement heating rate of 
2.4 times the flat-plate heating rate occurs for a gap width of 0.41 cm, whereas 
for the two smaller gap widths, the impingement heating was less than the sur- 
face heating. As the gap width increases, the flow in the longitudinal gap 
also influences the heating on the top surfaces. In addition, decreasing the 
slope of the forward-facing wall from 75® to 60® moderately decreases the 
impingement heating level and surface heating. 
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For turbulent boundary layers, the transverse flow is primarily three- 
dimensional in nature. The region affected by the three-dimensional flow 
extends at least three gap widths in the transverse and longitudinal directions, 

A maximum heating rate ratio of 4,5 occurs for a forward-facing wall slope of 
90^ and a gap width of 0.41 cm. The impingement heating generally increases 
with increases in gap width, gap length, and slope of the forward-facing tile 
wall. However, the effects of gap width and gap length are not completely under- 
stood because a combination of factors apparently creates a periodic fluctuation 
in the impingement heating for long gaps. In constrast, the geometric parame- 
ters have an overall opposite effect on heating of the surface of the tile imme- 
diately downstream of the impingement region; that is, if the parameter tends 
to increase impingement heating, the parameter tends to decrease the surface 
heating. The same periodic heating rate effect observed in the impingement 
region is evident on the top surface and occurred for both long and short gaps. 
Impingement heating rates can be reduced as much as 38 percent when the forward- 
facing wall slope is reduced from 90^ to 60^. In addition, the integrated heat 
load per unit width on the forward face is correspondingly decreased 23 percent 
even though the surface area is increased. 

Empirical relationships to correlate impingement heating with gap flow and 
geometric parameters for laminar and turbulent flow were developed through com- 
bined regression analysis and graphic techniques. These relationships accu- 
rately predict the impingement heating rate and the influence of gap flow and 
geometric parameters within the data range used to develop the relationships. 
Ninety percent of the impingement heating rates can be predicted within 
±29 percent. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
May 26, 1978 
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TABLE I.- TUNNEL CONDITIONS FOR CALIBRATION SERIES 


Run 

number 

^t,c» 

K 

Pt,c» 

MPa 

q, 

kPa 

R, 

per m 

deg 

Panel holder 
leading-edge 
configuration 

Boundary-layer 
flow conditions 

1 

1827 

4.2 

14.48 

1.15 X 10^ 

0.1 

Blunt 

Laminar 

2 

1861 

4.0 

13.99 

1.08 

4.5 

Blunt 

Laminar 

3 

1816 

11.9 

14.84 

3.19 

.1 

Blunt 

Laminar 








^Discrepancies In impingeinent region heating between runs 9 and 16 are not understood 
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TABLE IV.- INPUT FOR REGRESSION ANALYSIS 
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laminar; T, turbulent. 

y: XX ■ Run number, yyy = Thermocouple number 



Seven equally spaced at 6.4 — ► 5.1 

5.1 

-h + H- -h -I- + -I- — r 


+ + -f- -I- -H + + 

Five 

equally 

45.8 -1- + -h -h + spaced at 

I 8.9 



Figure 1.- Thermocouple locations on calibration panel. Dimensions are in cm. 
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Figure 2.- 60°/75° metallic tile array and geometric parameters. 





(b) 90° tile array. 

Figure 3.- Arrays of metallic tiles. 


L-78-104 





Figure 4.- Instrumentation for 60°/75° thin wall tiles. L = 15.24 cm; dimen 
sions are in cm. (For thermocouple designation where L = 30.48 cm, sub- 
tract 43 from each thermocouple number.) 




desi 








Dimensions are in 


Location I for 
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Figure 7.- Panel holder (with tile arrays installed) in Langley 
8-foot high-temperature structures tunnel. 






29 





Figure 9.- Typical history of events during model injection into test stream. 
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s> cm s, cm 


(a) Run 7: 0 = 60° and W = 0.18 cm. (b) Run 1: 6 • 75° and W = 0.18 cm. 



s, cm s, cm 


(c) Run 6: 0 = 60° and W = 0.41 cm. (d) Run 5: 0 = 75° and W = 0.41 cm. 

Figure 11.- Typical laminar convective heating rate distributions along 
longitudinal center line. Pt c “ ^*3 MPa; a = 0°; 6* = 0.57 cm; and 

L = 15.24 cm. ' 









s, cm 

(a) Run 9: 0 = 60^. 


s, cm 

(b) Run 16; 0 =* 75^, 



s, cm 


(c) Run 16: 0 = 90^. 

Figure 15.- Typical heating rate distributions in turbulent flow along 
longitudinal center line of tiles, ^ = 6.8 MPa; a = 7.5®; 

W = 0.18 cm; L = 15.24 cm. 










te contour plots in turbulent flow 
5°; L = 15.24 cm. 










0, deg 

Figure 19.- Effect of slope of forward-facing wall on impingement heating 
in turbulent flow. Pt, c “6.8 MPa; a - 7.5°; 6* = 1.01 or 1.36 cm; 

L = 1 5. 24 cm. 





W, cm 
(b) e = 90°. 


Figure 21.- Effects of 0, W, and L on surface he. 

(, = 6.8 MPa; a = 7.5°. 
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(a) Heating rate ratio plotted against 0. 
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Figure 23.- Comparison of empirical relationships and experimental data 

for turbulent flow. L = 15.24 cm. 
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W, cm 

(b) Heating rate ratio plotted against W. L = 15.24 cm. 


Figure 23.- Concluded. 


(o.Aipp) 


Laminar : 


Turbulent : 



/ ...I \ » 17 000 sin 6 

\ ^FP/ 


em 


/-i_\ = 5.58 (6*/W)-°*^^ (Ps/Patm)°*^^ ^ ® 

\^FP/„„ 


0.7C 


em 


Figure 24.- Empirical relationship for maximum impingement heating. 
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